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Abstract. We have obtained the first beryllium measurements of late F/early G-type stars in the old open cluster 
M67 (4.5 Gyr) and in the intermediate age cluster IC 4651 (1.7 Gyr). One member of the young cluster IC 2391 
(~ 50 Myr) was also observed. Our sample stars have effective temperatures within a range of -1-30 — 1-380 K from 
the solar temperature. All our sample stars, including the Sun and the young cluster star have, within the errors, 
the same Be abundance. This result implies that late F/early G-type stars undergo very little (if any) Be depletion 
during their main-sequence life-time. Since these stars have undergone some Li depletion, our finding is indicative 
of shallow mixing, i.e. of a mixing process that can transport surface material deep enough for Li burning to occur, 
but not deep enough for Be burning. As shown in previous studies, the Li vs. Be diagram is a powerful diagnostic 
of stellar interiors. In this context, we do not find any evidence of correlated Li and Be depletion; furthermore, 
a comparison with various models shows that the Be pattern of our sample stars is compatible only with models 
including gravity waves. This class of models, however, cannot reproduce the Li observations of M 67. 



Key words. Stars: abundances - Stars: Evolution 
Clusters and Associations: Individual: IC 4651 



Open Clusters and Associations: Individual: M 67 - Open 



1. Introduction 

The comparison between the chemical composition of me- 
teorites, assumed to be representative of the primordial so- 
lar cloud, and the Sun shows that lithium is 140 times un- 
derabundant in the solar photosphere (Anders & Grevesse 
1989). This result cannot be explained as due to the ac- 



tion of convection only since, according to standard stellar 
models, the convcctive zone of solar type stars on the main 
sequence (MS) does not reach deep enough layers to en- 
able lithium destruction (e.g., Pinsonneault 1997). Besides 
Li depletion in the Sun, several pieces of observational 
evidence have been found indicating that, in contrast to 
standard models, late F/early G-type stars do deplete Li 
while on the MS (e.g., Jeffries |2000[ Pasquini |2000| and 
references therein), witnessing the action of an additional 
process that can transport material from the surface to the 
stellar interior. Whereas several models have been devel- 
oped, which take into account a more complex physics, in- 
cluding mass loss, diffusion, slow mixing driven by rotation 
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or gravity waves (Swenson & Faulkner 1992, Chaboyer et 
al. 1995, Deliyannis & Pinsonneault 1997, Montalban & 
Schatzman 2000D , these models are still very poorly con- 
strained; in particular, the question remains how deep in 
the stellar interior the mixing extends. 

Since beryllium burns at temperatures 10^ K higher 
than Li (Li and Be are destroyed by proton capture at the 
temperatures of 2.5 and 3.5 xlO^ K, respectively) the si- 
multaneous determination of Li and Be abundances in the 
same star traces stellar mixing to different depths, thus 
providing stringent observational constraints on model 
predictions. Whereas until recently it was commonly ac- 
cepted that the Sun had a Be abundance a factor of 
about two lower than meteoritic, Balachandran & Bell 
( 1998| ) carried out a revised analysis of the solar spec- 
trum in the near-UV and found a solar Be abundance 
logn(Be)0 = 1.40, thus comparable with the meteoritic 
value logn(Be)= 1.42 (Anders & Grevesse [1989D . They 
concluded that Be is undepleted in the solar photosphere 
and suggested that mixing in the solar interior is shallower 
than previously thought. The result of Balachand ran fc 
Bell however is still controversial (e.g., King et al. 1997 ); 
in particular, the finding of little or no Be depletion in 
the Sun cannot be regarded as a general conclusion, since 
it relies on absolute abundance measurements, which are 
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Fig. 1. Comparison of the spectra in the Be (top panel) and Li (bottom panel) spectral regions of the stars S1252 
(solid line) and S1256 (dashed line) in the old cluster M 67. The figure clearly shows the lack of correlated Li and Be 
depletion; whereas a large difference in Li strengths is evident, the spectra in the Be region do not show a significant 
disagreement. Li equivalent widths of 33 and 11 mA were measured for S1252 and S1256 respectively, resulting in a 
factor of ~ 5 difference in Li abundance. 



hampered by large uncertainties. Furthermore, we do not 
know whether and to what extent our Sun is representa- 
tive of the whole population of stars with similar mass and 
evolutionary status. Finally, according to most of the pro- 
posed models, mixing depends on the whole evolutionary 
history of a star: hence, knowledge of the Be abundance 
for only one star at the age of the Sun does not allow 
discrimination between different models, since stars with 
different parameter evolution (e.g. rotation) may imdcrgo 
different amounts of mixing. Beryllium measurements in 



homogeneous samples of stars, such as stellar clusters, 
with well determined evolutionary status, chemical com- 
position, and age, arc thus crucial. In addition, since most 
models predict significant Be variations on long timescales 
(> 1 Gyr), a meaningful comparison between model pre- 
dictions and observations should include old stars. 

The Be feature that is best observed, the Be ii reso- 
nance doublet (A = 3130.420 A and A = 3131.064 A), lies 
close to the atmospheric UV cutoff; observations of Be arc 
thus very challenging and have so far been limited to the 
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star 


V 


(B-V)o 


Night 


Exp. 


S/N @ Be 


S/N @ 








(UT) 


(sec) 






a Cen A 


-0.1 


0.63 


2000 Apr 11 


4(b)/2(r) 


120 


340 


IC 2391 SHJM2 


10.3 


0.57 


2001 Feb 16 


2400 


30 


210 


M 67 S988 


13.18 


0.534 


2000 Apr 15 


6300 


40 


120 


M 67 S994 


13.18 


0.535 


2000 Apr 14 


10800 


50 


160 


M 67 S1252 


14.07 


0.587 


2000 Apr 16 


5400 


40 


110 


M 67 S1256 


13.67 


0.595 


2000 Apr 15 


9100 


45 


120 


M 67 S969 


14.18 


0.622 


2000 Apr 16 


5400 


40 


110 


IC 4651 EG7 


14.20 


0.557 


2000 Mar 26 


14400 


45 


180 


IC 4651 EG45 


14.20 


0.568 


2000 Mar 29 


10800 


40 


140 


IC 4651 AT2105 


14.02 


0.545 


2000 Mar 30 


9540 


40 


140 



Numbering for M 67 from Sanders (1Q77) 



Numbering for IC 4651 from EggenXroTl) and Anthony- Twarog et al. (1988) 





Fig. 2. The spectrum synthesis fit for the star S1252 in 
M 67. The histogram represents the observed spectrum; 
the dotted hne denotes the best fit, while the two dot- 
dashed lines indicate synthesis obtained with a factor of 
two lower/higher Be abundance. 



Fig. 3. Same as Fig. 2, but the spectrum synthesis for 
S1256 in M 67 is shown. 



while an abundance analysis including heavy elements will 
be presented in a forthcoming paper. 



brigh test stars observable in the field (e.g., Primas et al. 
2001 ) and in the close- by Hya des and Ursa Major clus- 
ters (Boesgaard & Budge 1989 , Garcia Lopez et al. 1995 , 
Boesgaard & King |00|). 

Thanks to the superior near- UV ca pabilities of UVES 
on VLT Kueyen (Dekker et al. 2000 ), we were able to 
acquire for the first time high resolution spectra of the Be 
region of 14*'' magnitude G dwarfs in old clusters. In the 
present paper we report the results concerning Li and Be, 



2. Observations 

Our sample includes five members of the solar age M 67 
cluster (4.5 Gyr), three members of the 1.7 Gyr old 
IC 4651, the solar-analog a Cen A, and one member of 
the young (50 Myr) cluster IC 2391. This star has just 
arrived on the zero age mai n seq uence and it is unde- 
pleted in Li (Randich et al. 2001 ); under the fully rea- 
sonable assumption that it is also undepleted in Be, it 
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Table 2. Effective temperatures, rotational velocities, and derived abundances 



star 


Tefl 


vsini 


EW(Li) 


log n(Li) 


log n(Be) 




(K) 


(km/sec) 


(mA) 






a Cen A 


5800 


2 


9 ±0.5 


1.24 ±0.11 


1.17± 0.11 


IC 2391 SHJM2 


5970 


9 


150 ±3 


3.45 ±0.09 


1.11 ± 0.13 


M 67 S988 


6153 


6 


12 ±2 


1.89 ±0.13 


0.88 ±0.13 


M 67 S994 


6151 


6 


25 ±2 


2.27 ±0.11 


1.16 ±0.13 


M 67 S1252 


5938 


4 


33 ±3 


2.35 ±0.11 


1.11 ± 0.13 


M 67 S1256 


5907 


4 


11 ± 1 


1.60 ±0.11 


1.16 ± 0.13 


M 67 S969 


5800 


4 


28 ±2 


2.06 ±0.10 


1.11 ±0.13 


IC 4651 EG7 


6061 


4 


52 ±3 


2.61 ±0.11 


1.11 ±0.13 


IC 4651 EG45 


6016 


4 


48 ±3 


2.52 ±0.11 


1.16 ±0.13 


IC 4651 AT2105 


6110 


8 


67 ±6 


2.82 ±0.12 


1.11 ±0.13 



Bell Bell 





3130 3131 
Wavelength 



3130 3131 
Wavelength 



Fig. 4. Same as Fig. 2, but the spectrum synthesis for 
AT2105 in IC 4651 is shown. 



Fig. 5. Same as Fig. 2, but the spectrum synthesis for 
the star in the young clusters IC 2391 is shown. 



can be used to set a zero point to our abundance scale 
and it may provide an indication of the initial Be abun- 
dance. No known binaries are included in our sample. The 
observations were carried out using UVES on VLT UT2 
(Kueyen). M 67 members and a Cen A were observed 
in Visitor Mode in April 2000; the young cluster mem- 
ber was also observed in Visitor Mode in Feb. 2001, as 
part of a different program; finally, the three targets in 
IC 4651 were instead observed in Service Mode during 
March 2000. All the stars were observed using the same 
set-up of UVES; namely, UVES was operated in Dichroic 
Mode using Cross Dispersers #1 and #3 in the Blue and 
Red arms, respectively. The Blue arm is equipped with 
a EEV 2048x4102 CCD, while the Red arm is equipped 



with a mosaic of two CCDs made by a EEV 2048x4102 
CCD and a MIT-LL 2048 x 4102 CCD. Such a combina- 
tion allowed us to cover the spectral ranges from ~ 3115 
to 3940 A in the blue and from - 4780 to 6810 A in the 
red. The 15/Ltm pixels and the use of a 1 arcsec wide slit 
(projecting into 4 pixels) and CCD binning 2x2 and 1x1 
in the blue and red, respectively, yielded resolving pow- 
ers R~ 40, 000 and R~ 45, 000. The spectra were reduced 
using MIDAS and following the usual steps: bias and inter- 
order background were subtracted from both science and 
flat-field frames; the science frames were then optimally 
extracted, flat-fielded, and wavelength calibrated. Finally, 
the extracted orders were merged. The exposure times of 
cluster stars ranged between 40 min and 3 hrs; the sample 
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stars were observed under different seeing conditions and 
variable airmasses, due to the long exposure times. The 
final S/N ratios per spectral bin are comprised between 
~ 30 and 50 in the Be region and between ^ 100 and 
200 in the Li region. A much higher S/N was achieved for 
a Cen A. The sample stars together with details on the 
observations are listed in Table 1. In Fig. 1 we show the 
comparison of the spectra of two M 67 members in the Be 
and Li spectral regions. 

3. Abundance and Error Analysis 

Effective temperatures were determined based on B-V col- 
ors using the same calib ration employed in similar studies 
(Soderblom et al. 1993 ). A surface gravity logg — 4.44 



and microturbulence ^ = 1.1 km/sec were assumed for 
all sample stars. Rotational velocities were retrieved from 
Pasquini et al. (in pr epara tion) for M 67 and IC 4651 and 
from Randich et al. (2001) for the young cluster member. 
Teff and vsini values are listed in Cols. 2 and 3 of Table 2. 

Li and Be abundances were derived in a consistent way 
using Kurucz model atmospheres and codes (ATLAS9, 
SYNTHE, and WIDTH9 - Kurucz |1993| ). Model atmo- 
spheres were computed for the exact stellar parameters of 
each star, by interpolating in the ATLAS9 grid of models 
officially released by Kurucz on CD-ROMs (Kurucz 1993| ) . 
Solar abundances were taken from Anders & Grevesse 

Be abundances were determined by spectrum synthe- 
sis using the line list extensively tested by Primas et al. 
( |1997D . This compilation has been constrained by using a 
large sample of stars, and after having run several tests 
on the possible presence of unidentified lines in the Be II 
doublet region. The main reason behind these tests (also 
carried out by similar studies) is the following: using the 
list of atomic and molecular lines derived from labora- 
tory measurements only and model atmospheres currently 
available, it is very difficult to obtain a satisfactory fit of 
the A3131.064 A line with the solar photospheric abun- 
dance (log ^(Be)© = 1.15, as determined by Chmielewski 



et al. 1975). Possible ways out include the introduction 
of one or more (still unidentified, hence "predicted") lines 
or adjusting the oscilla tor st rengths of neighbouring lines 
(cf. Garcia Lopez et al. 1995). Though different, these two 
approaches give similar results, a satisfactory fit of the so- 
lar spectrum with the abundance logn(Be)Q = 1.15. For 
our analysis, we adopted the line list from Primas et al. 
( |1997D , including a predicted Fe I line at 3131.043 A , the 
strength of which had been constrained by using a large 
sample of stars of different Teg, log 5, and spanning three 
orders of magnitude in metallicity (cf. the original work 
for more details). 

We started our analys is fro m fitting the Kurucz Solar 
Flux Atlas (Kurucz et al. 1984 ): using a solar model with 
Teff0 = 5777 K, log5o = 4.437, = 1.1 km/sec, we 
obtained a solar Be abundance logrt(Be)0 = 1.11, in very 
good a greem ent with the value found by Chmielewski 
et al. (1975). To gain confidence on the robustness of 



our analytical method, we then fitted our own spectrum 
of the solar-analog a Cen A (Toff= 5800 K, logg — 
4.44, [M/H]=+0.1, cf. Primas et al. |l997| ). We found 
log?i(Be)Q = 1.17, in good agreement with the Sun 
and with values derived from previous analyses (within 
the quoted uncertainties 
al. 



cf. Primas et al. 



1997, King 



et al. 1997). Then, we applied the same technique to 
all our sample stars, including also the member of the 
young cluster IC 2391. Examples of some spectrum syn- 
thesis fits are shown in Figs. 2-5, while all the Be (de- 
rived by demanding a satisfactory fit of hoth Be lines) 
and Li abundances are reported in Table 2, together 
with their respective uncertainties. For the records, us- 
ing model atmospheres computed with the overshooting 
option switched off (NOVER) translates into a systematic 
difference Alogn(Be) = —0.05 dex for all stars. 

Li abundances were derived from the measured equiv- 
alent widths of the Li i 6707.8 A feature (see Table 2). 
Given the resolution of our spectra and the rotational ve- 
locity of the sample stars, the Li feature was not blended 
with the nearby Fe i 6707.44 A line and, thus, we did not 
need to correct the measured Li EWs for the contribution 
of the latter. For completeness and consistency, Li abun- 
dances for stars from previous studies (namely, Jones et 
al. |1999| for M 67 and Randich et al. |2000| for IC 4651) 
have also been re-computed from the EWs published in 
the original works. We note that new abundances are in 
very good agreement with old ones. 

Errors in the derived abundances were computed from 
analyzing the dependence of Li and Be abundances to 
changes in the adopted stellar parameters. Uncertainties 
of ±100 K in Tcff , 0.2 dex in log g, and 0.2 km/sec in ^ were 
assumed. The final uncertainty listed in Table 2 include 
also the sensitivity of Be to the placement of the contin- 
uum (estimated to be on the order of ±3%) and of Li to 
errors in the measured EWs. There may be extra sources 
of systematic and random errors. In the case of Be, one 
may want to take into account extra uncertainties encoun- 
tered during spectrum synthesis. During our analysis we 
have established that including one "predicted" line (or 
adjusting the oscillator strengths of one or more lines) im- 
plies a systematic difference of ~ —0.05 dex for our sample 
stars here (as to say, without the predicted line, one has to 
increase the Be abundance by -1-0.05 dex in order to get a 
similar satisfactory fit). Another way is to run a synthesis 
including the entire Kurucz list of "predicted" (as opposed 
to "laboratory") lines, in order to quantify the global ef- 
fect on the continuum. Although related to the same prob- 
lem, this may be considered a random source of error. We 
found it affects the continuum level by ~ 3%, which cor- 
responds to ±0.05 dex in beryllium. Another weak point 
to be aware of is the well known discrepancy between the 
solar photospheric Be (logn(Be)0 = 1.15) and the abun- 
dance measured in meteorit es (lo gn(Be)0 = 1.42 deter- 
mined by Anders & Grevesse |l989|) . Though this may turn 
out to be onl y an apparent discrepancy (cf. Balachandran 
& Bell 1998 , who determined a solar photospheric Be in 
very good agreement with the meteoritic value by includ- 
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ing new continuum opacity sources), we decided not to en- 
ter the discussion which, as mentioned above, is stiU very 
much open, as also shown by AUende Prieto & Lambert 
(200C) who find that there is no missing UV opacity. 

In conclusion, based on the quality of our observations 
and on the fact that our main conclusions will be based 
on differential comparison between stars characterized by 
similar parameters, we believe that the relative error bars 
quoted in Table 2 are appropriate. In other words, thanks 
to our choice of targets of very similar nature and macro- 
scopic parameters, we can safely compare the relative light 
element abundances of our sample stars without being af- 
fected by systematic effects. This implies that differences 
in Be abundances are better constrained than absolute 
abundances. 



4. Results and discussion 

4.1. Be and Li abundances 

In Fig. 6 we show logn(Be) vs. Tcfr for our sample stars. 
The Sun is also included in the figure. Figure 6 and 
Table 2 clearly show that all our sample stars have, within 
the errors, the same Be abundances: in particular, we 
find for the Sun exactly the same abundance as for the 
50 Myr old IC 2391 member; stars in M 67 (with the 
possible exception of S988 that shows a somewhat lower 
abundance), IC 4651, and a Cen A also have a similar 
abundance as the Sun. We obtain an average abundance 
logri(Be)= 1.08±0.12forM67andlogn(Be) = 1.13±0.03 
for IC 4651. We stress that the finding of similar abun- 
dances is extremely robust since it is based both on direct 
comparison of the spectra and on a consistent analysis 
of stars with very similar macroscopic parameters. As we 
have mentioned, the Be abundance that we derive for the 
Sun is in agreement with th e valu e log n (Be) = 1.15 de- 
rived by Chmielewski et al. ( |l975|) ; this value, as well the 
abundances of our sample stars, are about a factor of two 
lower than the meteoritic abundance logn(Be) = 1.42. 
For the three clusters we have derived similar, about so- 
lar, metallicity (see also Randich et al. 2001 ): given the 
logn(Be) vs. [Fe/H] relationship (e.g., Boesgaard et al 
1999| ), this suggests a similar initial Be abundance. Two 



hypotheses can be proposed to explain the difference be- 
tween the abundances that we have derived for our sample 
stars and the meteoritic Be abundance: a) All our sample 
stars, including the young cluster star, have undergone a 
very similar amount of Be depletion; b) The difference be- 
tween our best fit Be abundances and the meteoritic value 
is simply due to different analytical methods adopted dur- 
ing the analyses. Under this hypothesis, none of our sam- 
ple stars has undergone any Be depletion, implying that 
stars in this temperature range including the Sun, do 
not burn Be while on the MS. Hypothesis a) would mean 
that some Be burning occurs during the pre-main sequence 
(PMS) phases, in contradic tion w ith all model predictions 
(e.g., Piau & Turck-Chieze 2001 ) and inconsistently with 
the observational evidence of no PMS Li depletion for stars 



in this Toff range (e.g., Randich et al. 2001 and references 
therein); in other words, the simultaneous Be depletion 
and Li non-depletion in the young star would be hardly 
explainable. For these reasons, we regard hypothesis a) as 
very unlikely. Hypothesis b) would instead be the proof of 
a shallow mixing process, i.e. a mechanism that is able to 
transport surface material down to the Li burning layer, 
but not deep enough to destroy Be. 

In Fig. 7 we plot logn(Li) vs. Tcff for our sample stars 
and th e star s studied by Jones et al. ( |1999|) and Randich 
et al. ( ^000 ). The present study confirms the star-to-star 



scatter in Li among M 67 members (Spite et al. 1987, 

Pasquini et al. |1997, Jones et 



1988 



Garcia Lopez et al 
al. 1999), indicating different amounts of Li depletion for 
stars with the same mass and chemical composition. In 
agreement with previous studies (Randich et al 



2000) 



stars in the younger IC 4651 cluster do not show any major 
dispersion, suggesting that the mechanism responsible for 
the dispersion observed in M 67 acts at relatively old ages. 



4.2. Be vs. Li depletion and comparison with the 
models 



As clearly pointed out by Deliyannis (200C), the 
logn(Be) vs. logn(Li) diagram is a very powerful tool to 
investigate the nature of stellar mixing and to discrim- 
inate between different models. It is worth emphasizing 
that in this diagram the shape of the relationship, i.e., 
the relative Li and Be abundances, rather than absolute 
abundances, is providing most of the information. In Fig. 8 
we show logn(Be) vs. logri(Li) for our sample stars and 
the Sun; the best fit relationship between logn(Be) and 
logn(Li) proposed by Boesgaard et al. (2001) for field 
stars with 6300 > Toff > 5850 K is also shown in the 
figure. 



Be surveys a mong field stars (e.g., Stephens et al. |1997 , 
Boesgaard et al. 2001 ) have reported the existence of a cor- 
relation between Li and Be depletion; this correlation was 
interpreted as the proof of the action of slow mixing driven 
by rotation; a similar result was found by Boesgaard & 
King (2002) for Hyades dwarfs warmer than 5850 K. At 
odds with the correlation found for field stars, and in con- 
trast with the prediction of slow mixing induced by rota- 
tion (see e.g., Deliyannis 2000| ), Fig. 8 does not show any 
clear evidence of correlated Li and Be depletion; stars that 
have suffered different amounts of Li destruction both in 
the same cluster and in different clusters, have, within the 
errors, the same Be content (see also Fig. 6). Our dat- 
apoints plotted in Fig. 8 do not follow the mean trend 
found for field stars with 6300 > Toff > 5850 K. This is 
not unexpected for the stars in our sample with Tcff be- 
low 5850 K (namely, the Sun, alpha Cen A, and S969 in 
M 67), for which no correlation between Li and Be is in- 
deed found for the Hyades. It is instead surprising that 
the warmest stars in our sample do not follow the same 
relationship as field stars. Most of the stars in M 67 and 
IC 4651 are only marginally in agreement with the rela- 
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6200 



6000 
Teff (K) 



5800 



6200 



6000 
Teff (K) 



5800 



Fig. 6. logn(Be) vs. effective temperature for our sample stars. Filled circles and triangles denote M 67 and IC 4651 
members, respectively, wfiile the open square indicates tlie young cluster member. The Sun (Sun symbol) and a Cen 
A (star symbol) are included in the figure. Curves denoting model predictions are shown in the figure. Dot-dashed 
curves indicate the predictions of models including gravity waves, with upper and lower curves denoting predictions 
at an age of 1.7 and 4.0 Gyr, respectively. Dotted and dashed curves indicate the predictions of rotationally induced 
mixing at 1.7 and 4 Gyr; for each curve style, upper curves refer to models with an initial rotational velocity of 10 
km/sec, while lower curves denote the predictions of models with an initial rotational velocity of 30 km/sec. In panel 
a) theoretical curves are normalized to an initial abundance equal to the meteoritic value logn(Be)= 1.42, while in 
panel b) the curves are normalized to the abundance logn(Be)= 1.11 that we derived for the young cluster member. 



tionship for field stars, with star S1256 in M 67 clearly 
showing much less Be depletion than expected; viceversa, 
the young cluster star has a Be abundance significantly 
below the expected value. If one changed the intercept of 
the curve to fit this star, all the other sample stars would 
then lie much above the relationship. We mention that 
different temperatures were reported in the literature for 
S1256, a couple of which are cooler than 5850 K. More 
in general, our te mper ature scale is different from that of 
Boesgaard et al. (|200l| ) and thus other stars in our sample, 
besides a Cen A, the Sun, and S969, might fall outside the 
interval 5850 < Tcff< 6300 K based on the calibration of 
Boesgaard et al. In order to further investigate this point, 
we re-determined the effective temperatures of field stars 
in the sample of Boesgaard et al. using our calibration. 
Whereas we do find differences as large as ~ 250 K, we do 
not find a systematic offset and the range 5850-6300 K of 
Boesgaard et al. reflects into an interval of 5810-6480 K 
in our scale. In other words, given the uncertainties in the 
temperature scales, the lower limit for the Li-Be correla- 



tion found for field stars cannot be defined in an absolute 
way. On the other hand, the discrepancy between the cor- 
relation found for field stars and the lack of correlation 
found for our sample stars cannot be explained as due to 
different temperature calibrations. 

In the following, our results are compared with the 
predictions of various models. 

Models including diffusion actually do not warrant a 
detailed discussion: as discussed by Deliyannis (2000) or 
by Stephens et al. (1997), diffusion implies simultaneous 
Li and Be depletion and thus does not fit the observed 
Be vs. Li diagram shown in Fig. 8. In Figs. 6 and 7, in- 
stead, we show a quantitative comparison with the pre- 
dictions of models including slow mixing processes; more 
specifically, we plot different curves denoting theoretical 
predictions by models including slow rotationally induced 
mixing (retrieved from Deliyannis & Pinsonneault 1997 ) 
and by models including gravity waves (kindly pr ovide d 
by J. Montalban -see also Montalban & Schatzman 2000). 
Following the two hypothesis a) and b) discussed in the 
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Fig. 7. logn(Li) vs. effective temperature for our sample stars. Symbols and curves are the same as in Fig. 6, but 
filled symbols refer to present data (namely, filled circles: M 67; filled triangles: IC 4651; filled square: young star; 
asterisk: a Cen A), while open circles and triangles denote Li abundances for M 67 and IC 4651 from previous studies 
(Jones et al. |1999| and Randich et al. |2000D . 
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denoted by the Sun symbol, a Cen A by a star symbol, and the young star in IC 2391 by an open square. Crossed 
symbols denote stars with Teff below 5850 K. The solid line indicate the best fit relationship found by Boesgaard et 
al. ( ^00 ID for field stars with TeS in the interval 5850-6300 K. 
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previous section, we show in panel a) of Fig. 6 model pre- 
dictions with an initial Be abundance equal to the mete- 
oritic value, as assumed in the original works presenting 
these models. In panel b) the curves are instead normal- 
ized to an initial abundance log?i(Be) = 1.11 as derived by 
us for the young, Li-undepleted star in our sample; under 
the assumption that this star has not depleted any Be, its 
abundance is indicative of the initial Be abundance in our 
abundance scale. 

Models including rotational mixing predict that some 
Be depletion has occurred already at 1.7 Gyr and that 
additional depletion should take place between 1.7 Gyr 
and the solar age. In addition, at a given age, they pre- 
dict a difference in both log7i(Li) and logn(Be) for stars 
with the same temperature, but different initial rotational 
velocities. The difference in Li abundance is much larger 
than the difference in Be abundance. In other words, ac- 
cording to these models, some dispersion in Be should be 
observed if cluster stars had different initial rotation rates. 
A larger dispersion in Li is predicted. 

The absolute Be abundances of M 67 and IC 4651 
members shown in panel a) are, within the errors, in fairly 
good agreement with the predictions by rotational mix- 
ing. However, as evidenced in Sect. 4.1, the young clus- 
ter member exhibits an abundance which is a factor ~ 2 
below meteoritic. The Be content of this star indeed ap- 
pears in agreement with the predicted Be abundance at 
an age between 1.7 and 4.5 Gyr, which is clearly incon- 
sistent with the much younger age of the star. Panel b) 
has been made under the assumption that the Be abun- 
dance of the young star would represent the initial Be 
abundance; therefore, and by definition, this panel does 
not evidence any inconsistency for the young cluster mem- 
ber. On the contrary, the observed Be abundances of both 
M 67 and IC 4651 members are not fitted by the mod- 
els including rotational mixing that predict larger deple- 
tion than observed. Only the warmest star in our sam- 
ple, S988 in M 67, is marginally in agreement with model 
predictions. Independently from the initial absolute Be 
abundance, rotational mixing is not able to reproduce 
the relative Li and Be abundances of our sample stars. 
A difference in Be abundance somewhat larger than our 
errors is predicted between 1.7 and 4.5 Gyr: we do not 
find evidences for such a difference. Furthermore, between 
6000 and 5900 K these models predict that a dispersion 
in logn(Li) of ^ 0.6 dex should correspond to a disper- 
sion in logn(Be) of ~ 0.2 dex, which is almost a factor 
of two larger than our errors. Stars S1252 and S1256 in 
M 67 show 0.75 dex difference in Li, but have virtually 
the same Be abundance, in contrast with the predictions. 
Even without considering the derived abundances, this re- 
sult is evidenced by the comparison of the spectra of the 
two stars in the Be and Li regions shown in Fig. 1. 

Gravity waves induced mixing predicts no significant 
Be depletion at 1.7 Gyr and only a small amount of de- 
pletion {~ 0.15 dex) for the coolest stars at the solar age. 
The predictions of gravity waves would be in contradic- 
tion with our results if hypothesis a) above is true and 



the sample stars have undergone some Be depletion; the 
curves showing the predictions of these models in panel a) 
indeed do not well fit the observed distribution. On the 
contrary, in panel b) models including waves provide the 
best fit to the observed Be vs. Teff diagram of our sam- 
ple stars. Mixing driven by waves also seems to reproduce 
fairly well the Li vs. temperature distribution of IC 4651 
stars, but is not able to reproduce the Li vs. tempera- 
ture morphology of the older M 67 shown in Fig. 7. In 
particular, available models cannot explain the dispersion 
in Li observed in M 67. Rotating models including grav- 
ity w aves are being developed (Montalban & Schatzman 
^OOOD ; whereas, different rotation rates could provide an 
explanation for the dispersion in Li, preliminary quali- 
tative analysis indicates that high rotation would result 
in lower predicted Li depletion thus not reproducing the 
lower envelope of M 67. 

We finally mention that we have not taken into account 
in this discussion Li depletion due to m ass los s. Models in- 
cluding mass loss (Swenson & Faulkner 1992 ) predict that 
Be depletion starts only when Li has been completely de- 
pleted and thus they would be qualitatively in agreement 
with our results. Mass loss, however, has been convinc- 
ingly shown not to work when applied to Li data alone 
(Swenson & Faulkner |l992|) . 



5. Conclusions 

We have obtained the first Be measurements of late 
F/early-G dwarfs in old open clusters. The Be abundances 
of old cluster stars are compared with the abundance of a 
young Li-undepleted star and with the solar abundance: 
we find that all the sample stars share, within the errors 
the same Be abundance, irrespective of their age and Li 
abundance. This in turn suggests that these stars, includ- 
ing the Sun, do not undergo any Be depletion during the 
main sequence, unless the very unlikely hypothesis is ac- 
cepted that they deplete Be, but do not deplete Li in the 
PMS phases. Our data are thus indicative of a shallow 
mixing process that extends deep enough in the stellar 
interior for Li burning to occur, but not deep enough to 
destroy Be. At variance with the results for field stars, 
we do not find any evidence of correlated Li and Be de- 
pletion, which would support shallow mixing driven by 
rotation. More in general, our observations and the com- 
parison with various models show that the understanding 
of the mechanism regulating the boundary layers between 
the convective zone and the radiative core is still not satis- 
factory and that the problem of depletion of light elements 
during the MS phases of late F/early G-type stars, includ- 
ing our own Sun, needs further investigation on theoreti- 
cal grounds. None of the models that we have discussed is 
able to simultaneously fit the observed Be vs. TeS and Li 
vs. Toff distributions of our sample stars and in particular 
their relative amounts of Li and Be depletion. Our sample 
stars cover a relatively small interval of effective tempera- 
tures; therefore we cannot extend our results/conclusions 
to warmer/cooler stars. Indeed, different mixing mecha- 
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nisms may be at work in stars much warmer/cooler than 
the Sun, with much thinner/thicker convective zones. In 
particiilar. we note that the possibility that the miying 



Pasquini, L., Randich, S., & Pallavicini, R., 1997, A&A 325, 
535 
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procesites at work in F a,nrl C-type stars are Hiftereiit 
is not unreahstic since the transition between G- and F- 
type stars corresponds to the region in the HR diagram 
where several changes occur; namely, the surface convec- 
tive zones gets thinner in F-type stars, possibly magnetic 
braking is not (or less) effective, and stars can retain high 
rotational velocities during their permanence on the main 
sequence. Most obviously, additional observations should 
carried out covering a larger temperature range. 
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